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RADIO FREQUENCY IMPEDANCE MAPPING 

Cross-Reference to Related Application 
This application clainxs liie benefit under 35 U.S.C. §119(e) of U.S. Provisional 
AppUcation 60/411,450, filed September 17, 2002, entitled ^TRADIO FREQUENCY 
IMPEDANCE MAPPING FOR MEDICAL IMAGESFG," by Sodickson et al, which is 
incorporated herein in its entirety. 

Field of the Invention 
The present invention relates generally to non-invasive imaging applications, for 
example, MRI. More particularly, the present invention relates to imaging techniques 
employing radio fi-equency (RF) coils to measure properties of a body bemg imaged. 

Background of the Invention 
A variety of imaging techniques have been employed to generate images of the internal 
characteristics of an object m diverse applications ranging fi^om medical imagmg to detection 
of prohibited materials in baggage at security checkpoints. For example, magnetic resonance 
imaging (MRI), X-ray computed tomography (CT), ultrasound, etc., have been widely utilized 
in medical applications to nnage the internal structures of an object, such as a patient. 
Generally, an imaging modality takes advantage of technology that facilitates discrhninating 
portions of an object based on one or more properties or characteristics of the object that can 
be measured. 

For example. X-ray CT includes measuring the attenuation of electromagnetic radiation 
as it passes througji an object. The resulting images carry information about the X-ray 
absorption characteristics of material withm the object, which is related to the atomic number 
of the material. MRI measures various magnetic properties such as relaxation times associated 
with various spin characteristics (e.g., realignment with an axis of precession) of material in a 
magnetic field. Ultrasound measures a material's capacity to reflect sound waves. 

While these techniques have been successful in discriminatmg between different 
materials (e.g., tissue, bone, etc.) to form images based on respective properties that can be 
isolated and measured, there may be subject matter of interest that cannot be readily 
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distinguished yfhen characterized according to the particular properties ejqjioited by a 
respective conventional imaging modality. For example, a tumor may have substantially the 
same X-ray absorption characteristics as the surrounding tissue, rendering the tumor effectively 
"invisible" to X-ray CT. 

Further, drawbacks of conventional imaging modalities such as X-ray CT and MRI 
mclude various practical limitations on flie equipment. MRI and X-ray CT devices are 
generally immobile and relatively expensive, often requirmg dedicated facilities for operation. 
The general bulk of such systems may prohibit the technology from being transported to a 
body, preventing these modaUties from being employed m the field or in emergency situations 
when a body (e.g., a patient) cannot be transported to the facility. 

Oflier imagmg techniques such as electrical unpedance tomography (BIT) have 
illustrated that dielectric properties (e.g., conductivity, permittivity, etc) of a body may be a 
viable diswiminating property to obtam images of internal structures of a body. In BIT, 
electrical currents are provided to the body through a set of electrodes applied to tiie surfece of 
the body, for example, the skin of a patient Change? in electric potential at each electrode are 
measured to determine dielectric properties of the body. However, BIT requires lhat electrical 
currents be applied dueotly to tiie body and that these currents be able to pass throughout the 
entire volume bemg imaged. Accordingly, some regions of a body may be difiBcult to image. 
For example, in human bodies, high resistance regions such as the skull may prevent brain 
images from being obtained. In addition, measurement uncertainties arise due to tiie 
impedance at the interfece between the electrodes and Ihe surface of the body being imaged. 

In magnetic mduction tomography (MIT), a pair of solenoid coils or gradiometer coils 
are positioned near an object to be imaged. The solenoid coils may operate as excitation coils 
and/or sensing coils. An excitation coil generates an oscillating magnetic field which is, in 
turn, detected by the sensing coil via magnetic mductance properties. The presence of a 
dielectric body between an excitation coil and a sensing coil perturbs the magnetic field sensed 
by the sensmg coil. Magnetic field perturbation results m a change in the mutual mductance 
between excitation and sensing coils (often measured as a phase shift in tiie magnetic field 
caused by eddy currents mduced in the dielectric object). Each excitation/sensing coil pair 
may provide a 'projection' of magnetic field perturbation. These projections may then be 
employed in various tomography techniques such as back-projection to reconstruct an image of 
the dielectric properties of the object. 
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However, MIT considers predominantly the properties of mutual inductance, for 
example, changes caused by eddy currents resulting from current induced in the dielectric 
object MIT utilizes solenoid coils or loop gradiometer pairs responsive predominantly to 
inductive coupling information. In general, loading effects on the inductive coupling between 
solenoid or gradiometer pairs are quite small and relatively difiBcult to measure. For example, 
a very precise measurement (on the order of a part in 1,000,000 or better) may be required to 
detect changes in inductive coupling resulting from the presence of a dielectric body. In 
addition, MTT relies on back-projection techniques which have been shown to yield 
quantitatively inaccurate results. 

Despite the shortcommgs of imaging modalities such as EIT and MIT, tiie dielectric 
properties of a body remain a viable characteristic for discriminating between regions of the 
body. For example, tumors generally have elevated values of both conductivity and 
permittivity relative to the surroundmg tissue. Hepatomas in rats, cancerous tissues in the 
breast, lung, colon, kidney and liver have been shown to have heightened dielectric properties 
relative to the tissue from which the cancer was derived. In addition, cerebral edema, 
spreading depression, myocardial ischemia and other pathologies may also exhibit dielectric 
contrast to the surrounding tissues and may benefit from imaging modalities capable of 
detecting this dielectric contrast (e.g., contrast in a material's conductivity and/or permittivity 
characteristics). 

As discussed above, MRI detects spin characteristics of target material to be imaged. 
MRI includes aligning the spin of nuclei of material bemg imaged in a generally homogeneous 
magnetic field and perturbing the magnetic field with periodic radio frequency (RF) pulses in 
order to measure the nuclear magnetic resonance (NMR) phenomenon of the material being 
imaged. To invoke the NMR phenomenon, one or more resonant coils are provided that 
genemte the RF pulses at a resonant frequency that matches a Larmor frequency (i.e., the rate 
at which a nucleus precesses about an axis) of certsun tissue in order to excite the nuclei such 
that Ihey precess about an axis in the direction of the applied RF pulse. When the RF pulse 
subsides, the nuclei realign with the magnetic field, releasing energy that can be measured. 

However, when a resonant coil is placed in proxunity of a load, for example, a patient 
or other object to be imaged, various properties of the resonant coil may be affected. In MRI, 
this loading effect tends to negatively impact the operation of the device by altering the 
resonant frequency of the coil and causing other generally undesirable changes in the in coil 
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properties. This loading effect depends in part on the dielectric properties of the load. 
Changes in resonant frequency of the coil may reduce the device's ability to excite the nuclei 
of the material beuig imaged (e.g., by creating a mismatch between the coiPs resonant 
frequency and the Larmor frequency of the target material) and negatively impact the quality 
of the resulting images. Tlie effects of coil loading complicate MRI to the extent that resonant 
coils are often tuned or adjusted to compensate for the generally undesirable loading effect 
caused by the body being imaged. 

Summary of the Invention 

One embodiment according to the present invention includes a method of determining 
dielectric properties of a body positioned proximate an array of coils having one or more 
resonant properties, the method comprising acts of detecting a change in at least one resonant 
property of at least one of the coils in the array, determining at least one electromagnetic 
property of at least one region of the body from the change in the at least one resonant property. 

Another embodiment according to the present invention includes a method of 
detenuitiing dielectric properties of a body, the method comprising acts of positioning the body 
proximate a plurality of coils, measuring at least one property of at least one of the pliirality of 
coils, and determining at least one electromagnetic property of at least one region of the body 
from the at least one property based on at least two of a resistive coupling, a capacitive 
coupling, and an inductive coupling between at least two of the plurality of coils. 

Another embodiment according to the present invention includes an apparatus for 
determining dielectric properties of a body. The apparatus comprises a plurality of coils having 
one or more resonant properties, a first component coupled to the plurality of coils and adapted 
to provide at least one measurement of the plurality of coils indicative of a change in at least one 
resonant property of at least one of the plurality of coils, and a second component coupled to the 
first component to receive the at least one measurement, the second component adapted to 
determine at least one electromagnetic property of at least one region of the body based on the 
change in the at least one resonant property. 

Another embodiment according to the present invention include an apparatus for 
determining dielectric properties of a body. The apparatus comprises a plurality of coils, a first 
component coupled to the plurality of coils, the first component adapted to provide at least one 
measurement of at least one property of the plurality of coils, and a second component coupled 
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to the first component to receive the at least one measurement, the second component adapted to 
determine at least one electromagnetic property of at least one region of the body from the at 
least one measurement based on at least two of a resistive coupling, a capacitive coupling, and 
an inductive coupling between two or more of the plurality of coils. 

Another embodiment according to the present invention includes a computer readable 
medium encoded with instructions capable of being executed on at least one processor, the 
instructions, when executed by the at least one processor, performing acts comprising defining 
an electromagnetic model of the coil array, receiving an input including a measured impedance 
matrix of the coil array, logically partitionmg a volume associated with the model of the coil 
array and the body into a plurality of regions, assigning trial values respectively to each of the 
plurality of regions, the trial values includuxg at least one of conductivity, permittivity and 
permeability, generating a trial impedance matrix from the assigned trial values according to the 
electromagnetic model of the coil array, and reducing a distance between the trial impedance 
matrix and the measured impedance matrix. 

Brief Description of the Drawings 
FIG. 1 illustrates a loading effect by mutual inductance between a conducting loop and 
a dielectric loop; 

FIG. 2 illustrates one embodiment according to the present invention of a resonant coil 
suitable for arrangmg in an array to determme dielectric properties of a body to be imaged; 

FIG. 3 illustrates one embodiment according to the present invention of a method for 
determining dielectric properties of a body by measuring the loading effect of the body on a 
coil array; 

FIGS. 4A and FIGS. 4B illustrate how spatial mformation on the dielectric distribution 
of a volume may be obtained using couplings of magnetic flux between coils in an array; 

FIG. 5 illustrates one embodiment according to the present invention of a method for 
determining dielectric properties of a body being imaged by measuring an impedance matrix of 
a coil array; 

FIG. 6 illustrates a coil array and a body that may be employed in connection with the 
method illustrated in FIG. 5; 

FIG. 7 illustrates one embodiment according to the present invention for calibrating a 
coil array; 
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FIG. 8 illustrates one embodiment according to the present invention of a coil array and 
a schematic of a model of the array; 

FIG. 9 illustrates one embodiment according to the present invention of an image 
acquisition system adapted to obtain one or more RFIM images; 

FIGS. lOA, lOB and IOC illustrate several embodunents of matching circuitry Ihat may 
be used to measure properties of a coil array according to various aspects of the present 
invention; 

FIG. 1 1 illustrates one embodiment according to the present invention of a method of 
acquiring one or more RFIM hnages; 

FIGS. 12A and FIGS. 12B illustrate one embodiment according to the present invention 
of an arrangement of a coil array; and 

FIG, 13 illustrates another embodiment according to the present invention of an 
arrangement of a coil array. 

Detailed Description 

While loading effects may cause generally undesimble results in conventional imaging 
modalities such as MRI, Applicant has recognized that the sensitivity of resonant coils to 
loading effects may be exploited to effectively measure dielectric properties of a body 
proximate to the coil. In particular. Applicant has identified and appreciated that a change in 
resonant properties of one or more resonant coils due to loading may provide information 
about the distribution of dielectric characteristics of the loading body. 

The term "resonant property" refers to a characteristic, trait or feature mdicative of the 
resonance of a coil including, but not limited to, resonant frequency, phase, decay 
characteristics, lineshape, etc. The term "loading effect" refers generally to any change in the 
electrical or magnetic property of an resonant coil resulting from the presence of a body. For 
example, the loading effect may mclude, but is not limited to, changes in the resonant 
properties of the coil, electromagnetic fields generated by the coil, coil impedance, etc. 

In one embodiment according to the present invention, changes in resonant properties 
of one or more resonant coils may be measured to determine the dielectric properties of a body. 
For example, impedance characteristics of one or more coils resulting from shifts in the 
resonant frequency of the coils may be measured to determine the distribution of dielectric 
characteristics of a loading body. Dielectric properties or characteristics refer generally to 
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electrical conductivity, electrical permittivity and/or magnetic permeability (also referred to 
simply as conductivity, permittivity and permeability, respectively). The term "bod/' will be 
used herein interchangeably with the term "objecf to refer generally to any mass capable of 
modifying one or more properties of a resonant coil, for example, by acting as a load. 

A resonant coil typically resonates at frequencies within the radio frequency (RF) 
spectrum. The RF spectrum is typically considered to span frequencies from approxunately 3 
kHz to lOOGHz. However, a resonant coil may be made to resonate outside this range. A 
resonant coil may be characterized by an inductive-capacitive-resistive (LCR) circuit. For 
example, a resonant coil may be comprised of one or more series capacitors connected by 
conducting material to form a loop having an inductance L, a resistance R, and a capacitance 
C. Accordingly, the impedance of such a coil at a frequency m may be expressed as, 



When such a coil is placed near a loading body, the impedance characteristics of the 
coil may be affected. Changes in impedance characteristics may be measured to obtain 
information about the dielectric make-up of the body loading the coil. For example, a body 
acting as a load on the coil may be characterized as one or more dielectric loops as illustrated 
in FIG. 1. Dielectric body 1 10 may also be considered as an LCR circuit having an inductance 
/, capacitance c and resistance r. The values of /, c and r may depend on the dielectric 
properties of the body such as conductivity and permittivity. 

When dielectric body 110 is placed proximate to a resonant coil 100 and the coil is 
operated, the resonant electromagnetic (EM) field generated by the coil 100 induces a current 
in dielectric loop 110, which in turn effects the properties of the resonant coil. For example, 
when a voltage V is provided to coil 100 such that an applied current 102 is generated in coil 
100, dielectric body 110 experiences an induced current 104. Accordingly, the resonant coil 
and dielectric loop inductively couple, that is, they exhibit a mutual inductance. In the 
presence of a dielectric body, the impedance of the RF coil at frequency co may be expressed 
as. 



(1) 
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That is, the impedance of the coil is affected by the presence of the load and depends in 
part on the dielectric properties of the loading body. Equation 2 provides one description of 
the impedance of a resonant coil in terms of the properties of a dielectric body. 

Applicant has identified and developed methods of determining dielectric properties of 
a body jfrom the impedance characteristics of one or more resonant coils operated in the 
presence of tiie body. The term "operate" with respect to a coil, refers generally to directiy 
generating a current (e.g., from a power source) in the conductive loop of the coil. For 
example, a voltage or current source may be provided at an input to operate the coil. The term 
"radio frequency impedance mappmg*' (RFIM) will be used herein to describe generally any of 
various methods of mapping impedance characteristics of one or more resonant coils to 
dielectric properties of a body coupled to the one or more resonant coils. 

When a resonant property of a coil changes, for example, due to the presence of a 
loading body, the impedance characteristics of the coil experience a corresponding change. In 
some exemplary coils, a resistance may change from .S Ohms on resonance when the coils are 
not loaded to 3 Ohms or more when loaded. In addition, eddy currents and displacement 
currents may modify the inductance and/or capacitance of such coils by an amount sufGcient to 
shift a resonant frequency of a respective coil by one percent or more. 

Accordingly, changes in resonant properties of an array of coils due to the presence of a 
loading body may be determined by measuring the impedance characteristics of the array to 
obtain information about the dielectric distribution of the body. When a coil is operated at its 
resonant frequency, an impedance of the coil is quite small and loadmg elffects may be quite 
large by comparison. As a result, changes in resonant properties may facilitate measurement 
precision that is far less demanding than other techniques, for example, MIT. 

The term "impedance characteristics" refers generally to measurements from which 
impedance may be derived. For example, impedance characteristics may be direct impedance 
measurements, resistance values, or may be measurements of scattermg parameters (S- 
parameters) which can be converted to and from impedance values. The term "array" refers 
generally to a plurality of coils arranged in a generally known relationship to one another. 

It should be appreciated from the foregomg that determining dielectric properties of a 
body may include measuring one or more impedance characteristics of a coil array. FIG. 2A 
illustrates one embodunent according to the present invention of a resonant coil that may be 
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arranged in an array and employed to obtain dielectric properties of a loading body. Coil 200 
may be a substantially planar coil capable of radiating electromagnetic energy. Coil 200 may 
include a generally non-conductive substrate 210 on which to mount a conducting loop. The 
conducting loop may comprise a plurality of conductive strips 220a-220d coupled together by 
capacitors 230a-230d. Hie conductive strips and the capacitors together form a resonant LCR 
circuit. The conductive strips may be formed o^ for example, copper material or any other 
conductive metal or substance. 

Coil 200 may be constructed, for example, from a printed circuit board ^PCB) board 
having a layer of copper over a substrate of a plastic material or other insulative material. The 
copper may then be etched to form the conductive loop. Coil 200 may also be constructed by 
applying copper strips to an adhesive substrate, or by any other method suitable for providing a 
conductive loop. A coil constructed of multiple turns or layers of conductor may be used. 
Capacitors 230a-230d may be any of a variety or combination of lumped elements or, in some 
embodiments, may be suitable distributed elements. Numerous combinations and variations in 
implementing a conductive loop will occur to those skilled in the art. The variety of methods 
of implementing a resonant circuit is considered to be within the scope of the invention. 

One of the capacitors (e.g., capacitor 230d) may be a matching capacitor. The 
matching capacitor may be used to match the input impedance of the coil to an RF power 
source. Matching capacitor 230d may be placed between leads of a connector adapted to 
couple to a power source capable of generating a current in the conductive loop of the coil. 
For example, coil 200 may be connected to a voltage source capable of providing a broadband 
electric field pulse via a Bayonet Neill Concelman (BNC) cable or other suitable coaxial cable 
240. Matching capacitor 230d may be chosen so that the coil has an input impedance that 
matches the impedance of the cable (e.g., a 50 Ohm or 75 Ohm cable). It may be desirable to 
match the input impedance of the coil at the nominal resonant frequency of the coil. It may 
also be desirable to match the input impedance of the coil at its resonant frequency when it is 
loaded by a body having similar characteristics to the body to be imaged. 

The term '^nominal" in connection with a resonant frequency refers to the frequency at 
which liie coil was adapted to resonate. In particular, the nominal frequency refers to the 
resonant frequency of an unloaded coil. As such, a change in resonant frequency refers to the 
change from the nominal value of the resonant frequency of the coil. Likewise, a shift in 
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resonant frequency refers to a change from a nominal resonant frequency that may occur, for 
example, when a resonant coil is loaded. 

The nominal resonant frequency of coil 200 depends in part on the arrangement and 
values of capacitors 230a-230d. FIG. 2B Ulusttates a circuit schematic of one embodiment of a 
resonant coU according to flie present mvention. Circuit 201 may schematically represent a 
coil similar to coil 200 having capacitance values as indicated. Specifically, capacitors 230a'- 
230c* may be 30 pF capacitors, and the matching capacitor 230d' may be a 94 pF capacitor. A 
coil arranged substantiaUy as shown in FIG. 2B may have a resonant frequency of 
approximately 110 MHz. It should be ^reciated that coils configured to have different 
nominal resonant frequencies may have differently valued capacitances. Similarly, coils 
having a nominal resonant frequency of 1 10 MHz may be designed using components different 
from the exemplary coil described above. 

Various properties of the loading effect may be ejqjloited using resonant coils that 
resonate at any number of fi^quencies. At higher resonant fi«quencies, for example, at 
frequencies approachmg I GHz, the loading effect is relatively strong. However, the 
electromagnetic fields generated at such high frequencies may have difficulty penetrating 
material comprising the body to be unaged. For example, when the body to be imaged is a 
human patient, high fi«quency EM fields may have difficulty penetrating tissue due in part to 
the so-called "skin depth" effect. 

At relatively low fitsquencies, for example, frequsmcies well below 100 MHz, the 
contrast between dielectric properties of a body may be enhanced. However, at relatively low 
frequencies the loading effect tends to becomes less significant and more difiGcult to measure. 
However, any resonant frequency may be used and may be chosen m consideration of Ae 
properties of the body being imaged and the characto-istics of a given coil or coil array. 

In one embodiment, coU capacitors (e.g., capacitors 230a-230d) may be replaced with 
varactors having adjustable capacitance. Accordingly, a coil may be provided having a 
variable resonant frequency. Such a coil array may be adjusted to resonate at one or more 
different frequencies in consideration of the type of body being imaged, the arrangement of 
coils, or other considerations of the imaging enviromnent. 

It should be appreciated that HGS. 2A and 2B illustrate only one of many varieties of 
resonant coils that may be suitable for practicing various aspects of the present mvention. 
Other coils may be suitable, includmg, but not limited to, birdcage coils, surface coils, volume 
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coils, transmission line array coils, various other resonant coils, etc. A resonant coil may 
comprise any number and arrangement of reactive components (e.g., capacitors and inductors) 
by which resonant electromagnetic fields may be generated. 

For simplicity and clarity, various aspects of the present invention will be illustrated in 
connection with the generally planar coil illustrated in FIG. 2. For example, FIGS. 6 and 9 
illustrate coil arrays comprising a plurality of coils schematically depicted as the general 
conductive loop connected by lumped capacitors. However, this depiction is meant to indicate 
the presence of a resonant coil and not the specific type or arrangement (e.g., planar, birdcage, 
etc.). As such, any type of resonant coil could be used. The invention is not limited to the 
coils specifically illustrated herein and contemplates use with any of the coils mentioned above 
or any other coils capable of generating resonant EM fields. 

FIG. 3 illustrates one embodiment according to the present invention of a method for 
determining dielectric properties of a body acting as a load on an array of RF coils. The coil 
array may comprise a plurality of coils each having a resonant frequency. That is, an unloaded 
coil may be configured to resonate at a generally known nominal frequency. Each coil in the 
array may be designed to have a same or different nominal resonant frequency. 

In step 3 1 0, a body (step 305) to be imaged may be positioned proximate to the 
plurality of coils m the array. In particular, the body is placed in a spatial relationship with the 
array such that the body acts as a load, causing some measurable loading effect on the array. 
The body is typically positioned in a generally known relationship with the array. The body 
may comprise one or more regions with unknown dielectric properties. For example, the body 
may be a patient having regions of homogeneous and inhomogeneous regions of conductivity 
and/or pemiittivity, or the body may be baggage such as airlme luggage to be analyzed for 
items having certain conductivity and/or permittivity characteristics. 

In step 320, a change in one or more resonant properties the coil array may be detected. 
For example, one or more measurements of properties of the coil array may be obtained. As 
described above, the loading effect resulting from placing the load proximate to the one or 
more coils will affect various properties of the coils. For example, various impedance 
characteristics of the coil array may be measured by operating one or more of the coils and 
measuring the impedance characteristics of one or more other coils in the array. 

In step 330, the change in the one or more resonant properties of at least one of the 
coils may be used to deteraiine dielectric properties of the body. For example, one or more 
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properties measured from the coil array may be employed to determine the dielectric properties 
of the body. Measured impedance characteristics may indicate properties of the load, for 
example, as expressed m the equation 1. Impedance characteristics of a loaded coil array may 
encode information about resonant frequency shifts, and resistive, capacitive, and inductive 
coupling between coils m the array that can be used to determine the dielectric distribution of 
the load. 

It may be desirable to form one or more images of Ihe dielectric properties obtained 
from a body. An image typically represents intensity as a function of space. The term 
"intensity" refers generally to a magnitude, degree and/or value at some location m the image, 
which in turn corresponds to a respective region of space. 

For example, in an X-ray image, intensity generally represents the absorption 
characteristics (e.g., density) of scanned material at a particular location in space. An image 
(e.g., an RFIM image) obtained via methods accordmg to the present invention may have 
intensity values representing conductivity, permittivity, permeability and/or other 
electromagnetic characteristics (e.g., electric and/or magnetic field magnitude, direction, phase, 
or some combination of the above) as a function of space or may have an intensity representing 
some combination of one or more of the above. In a two dimensional image, each intensity 
value may correspond to an area in space and is referred to as a pixel. In a three dimensional 
image, each intensity value may correspond to a volume in space and is referred to as a voxel. 

Accordingly, an RFIM image may involve determining spatial information about the 
dielectric properties of a body, that is, a dielectric distribution of the body may be obtained 
from impedance measurements of the one or more resonant coils. In one embodiment 
according to the present invention, an array of resonant coils is employed to encode spatial 
information about the dielectric properties of a body. 

In another embodiment, a single coil may be used to determine dielectric properties of 
the body. The single coil may be translated, rotated or otherwise varied in location with 
respect to the body being imaged. One or more properties of the single resonant coil may be 
measured at each of the locations. Accordingly, each placement of the coil will provide 
information about the overlap of the coil's electromagnetic field with the dielectric variation 
inside the body to be imaged as described in further detail below. 

FIGS. 4A and 4B illustrate principles by which an array of coils may be employed to 
obtam spatial mformation about the dielectric properties of a body. Array 450 includes a 
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plurality of coils 400a-400L In FIG. 4A, some of the magnetic field lines generated fi-om coil 
400a when the coil is operated are shoWn as magnetic field lines 405b-405L The partial 
magnetic field illustrated is merely schematic and represents a subset of the magnetic field 
lines that may be generated by coil 400a. 

The exemplary magnetic field lines shown are employed to illustrate that, when coil 
400a is operated, each of the other coils may experience a non-zero magnetic flux as a result of 
the magnetic field generated by coil 400a. The magnetic flux through each coil may induce a 
current in the respective coil proportional to the magnitude (i.e., the rate of change) of the flux. 
Accordingly, this induced current may be measured. While only magnetic field lines for coil 
400a are illustrated, it should be appreciated that each of the coils, when operated, will 
generate a magnetic field that will induce a current in other coils in proportion to the magnetic 
flux passing through the coil. 

In FIG. 4B, a dielectric body 410 is positioned within the array. The dielectric 
properties of the body 410 may affect the magnetic flux passmg through one or more of the 
coils as a result of the operation of coil 400a. For example, the conductivity of body 410 in the 
presence of the magnetic field generated by coil 400a may induce eddy currents that oppose 
the magnetic field applied by coil 400a, affecting the magnitude of the flux through one or 
more of the coils. The change in magnetic flux, in turn, affects the current induced in the 
respective coil. 

For example, magnetic field lines 405h and 405i may be substantially altered by the 
presence of body 410. Similarly, proximate magnetic field lines 405g and 405j may also be 
affected by the dielectric properties of body 410. Other magnetic field lines generated by coil 
400a may be less or even negligibly affected by the presence of body 410. Accordingly, which 
coils experience a change m magnetic flux and the extent of the change provides an indication 
of the location of tiie dielectric body 410 and its dielectric properties. 

For each operating coil, a pattern of magnetic flux perturbation (which results in a 
modification in induced current) will resuft m each of the other coils as a function of the 
location and dielectric properties of body 410, Since the location of each coil and the baseline 
current response of each coil in the array may be known, the patterns of flux modification 
yields information about the spatial distribution of a body's dielectric properties. 

Applicant has identified and appreciated that in addition to the magnetic mduction 
component (i.e., the inductive coupling between coils), the electric fields of the coil array 
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provide another mechanism for determinmg the spatial distribution of dielectric properties of a 
body. The electric field of a coil results in part from potential differences across inductors and 
capacitors and also jfrom any time-varying magnetic fields of the coil. These electric fields 
cause a current in conducting bodies and displacement currents in insulating (but polarizable) 
bodies. Displacement currents may modify the reactance of the coil, while currents in 
conducting bodies may modify coil resistance. Accordingly, the interaction of a coil array also 
has an electric field component (e.g., resistive and capacitive coupling between coils) that may 
be exploited to determine the dielectric distribution of a body. 

In MIT, a back-projection algorithm similar to techniques widely used in X-ray CT is 
employed to determme dielectric characteristics of a body. For example, a solenoid coil and a 
gradiometer may be arranged such that a body to be imaged is disposed between them. A 
mutual inductance between the solenoid coil and the gradiometer may be measured and the 
solenoid/gradiometer pair then shifted relative to the body and another measurement taken. 
This process may be repeated, for example, around a circumference of the body to provide a 
series of "projections" of inductive coupling between tiie solenoid/gradiometer pair. However, 
back-projection has been shown, in some cases, to yield quantitatively and qualitatively 
inaccurate results. In MTT, these errors may be exacerbated by relatively poor sensitivity to 
loading changes of non-resonant detectors. 

While back-projection may also be used to in the fi-amework of resonant coils 
according to various aspects of the present invention, it may have lunitations in extractmg 
spatial information about the distribution of dielectric properties of a body. Applicant has 
recognized and developed methods of determinmg the dielectric distribution of a body that 
incorporates the magnetic component and the electric component of an array of resonant coils. 
In particular. Applicant has developed methods incorporating effects of resistive coupling, 
capacitive coupling and inductive coupling between resonant coils in an array. 

Applicant has recognized that disturbances in the electric and magnetic fields caused by 
loadmg a coil array may be quantified using Maxwell's equations. In particular. Maxwell's 
equations may be employed to compute an impedance matrix of the coil array. The term 
"impedance matrix" refers generally to any ordered, related or otherwise correlated set of 
impedance characteristics of one or more coils in the presence of one or more other coils 
and/or dielectric bodies. An impedance matrbc may include any impedance characteristics, for 
example, S-parameters or other measurements fi-om which impedance may be derived. 
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A dielectric body may have spatially varying electric conductivity, permittivity and/or 
permeability. For example, a body may have an electrical conductivity o-(3c), electrical 
permittivity (3c) , and magnetic permeability //f3c; where the vector x is a direction of the one or 
more axes over which the conductivity and permittivity vary. Variations in magnetic 
permeability p, may be small in biological tissues, in which case, the permeability may be 
represented by a constant value fio. However, some materials, such as metals, may have 
appreciable and measurable variation in magnetic permeability. 

It should be appreciated that electromagnetic field energy is reduced when work is done 
on matter within and under the influence of the field. For example, when current is induced in 
a dielectric body, the energy in the inducing electromagnetic field is reduced m an amount 
proportional to the energy of the induced current. Accordingly, conservation of energy 
principles may be employed to ascertain various characteristics of a dielectric body by 
measuring the energy in the EM fields "lost" to the dielectric body. 

For example, consider an array of N coils arranged m a substantially known 
relationship to one another. Let E,{x) and b,{x) denote, respectively, electric and magnetic 
fields resulting fi-om a unit current m the coil of the array. If the currents and fields are 
time-varying (e.g., of the form c"'*') an impedance matrix for the array may be expressed as, 

2(^=j|^W^/(4^/(^)-'-«[^W^;(^>ty( (3) 

where Zy is the impedance of the coil in response to operating the /' coil. The first 
term of equation 3 gives the resistive loading of the array fi-om ohmic losses and the second 
and third terms give the capacitive and inductive loading of the array, which modify the array's 
reactance. The last term incorporates the effects of energy radiated out of a volume of interest. 
Stated differently, equation 3 describes loading effects of a body by considering each of 
resistive, capacitive and inductive coupling between the coils in the array, expressed in terms 
of coil impedance. 

The first integral in equation 3 may be evaluated for some suitably large volume 
surrounding the coil array. Radiation losses, which are represented by the final term in 
equation three, may be incorporated by integrating the Poynting vector (i.e., the cross product 
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of the electric and magnetic fields) over a surface that describes a volume of interest. This 
volume of interest, referred to as tixe "imaging volume", typically includes at least a portion of 
the loading body. The imaging volume, therefore, refers generally to the portion of space 
where it is desired to ascertain the dielectric distribution, for example, the distribution of 
conductivity and/or permittivity characteristics. 

Equation 3 illustrates some of the spatial encodmg mechanisms described above. For 
example, the first term of equation 3 (i.e., the resistance matrix) measures the overlap of the 
body's conductivity with a product of the electric fields of the i**" and coil as a function of 
space. Spatial encoding by mutual inductance is implicit in equation 3. That is, the electrical 
and magnetic fields themselves depend on the spatial distribution of the dielectric properties of 
the body. For example, the "blocking" or "screening" of the magnetic field by eddy currents as . 
described in connection with FIGS. 4A and 4B illustrate one way in which the electromagnetic 
fields depend on tiie dielectric properties of the body being imaged. 

It should be appreciated that inductive coupling between coils is fully expressed in the 
final term. This is in contrast to the linear approximation and back-projection techniques used 
in MIT. In fiirttier contrast, equation 3 fully expresses the resistive and capacitive coupling 
between the coils, providing additional information fi-om which determinations about the 
dielectric distribution within an imaging volume may be obtained. By considering both the 
electric and magnetic field components of a coil array, various aspects of the present invention 
may provide a fiill description of the relationships between coils in an array in the presence of 
a loading body. 

Equation 3 expresses the impedance of a coil array in terms of the electromagnetic 
fields generated by the array and the dielectric distribution of the imaging volume (e.g. a 
volume including a loading body>based on conservation of energy principles. However, it 
should be appreciated that there may be other ways of expressing coil array impedance 
characteristics in terais of the dielectric properties of a body that are suitable for implementing 
various aspects of the present invention and other formulations may be used witiiout departing 
from the scope of the invention. 

For example, dielectric properties of a body may be determined by computing currents 
that result from applying a known set of voltages to an array of resonant coil (e.g., using a 
finite difference time domain (FDTD) algorithm as described in more detail below). For 
example, a matrix of induced current values Iij(t) may be computed, where i indicates the coil 
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where a known voltage V is applied and j indicates the coil at which an induced current is 
being computed. 

By using Ohm*s law and by taking the Fourier transform of the current and voltage 
data, the intrinsic* impedance of the coils may be expressed as. 



where Fis a diagonal matrix containing the known voltages. The input impedance of 
the coil array may then be computed by combining Zintrinsic with the impedance of the matching 
circuitry. For example, in coils similarly arranged as shown in FIG. 2, the parallel circuit 
computation may be expressed in matrix form as. 



However, other forms may be used as well and may depend on the arrangement and 
topology of the coils. 

It should be appreciated that by solving Maxwell's equations within an imaging 
volume, not only may dielectric properties of a loading body be determined but other 
electromagnetic properties as well. For example, properties of the electric and/or magnetic 
fields may be determined eis a function of space. Accordingly, an image may be formed having 
intensity values indicating a magnitude, a direction, and/or phase of the electromagnetic 
environment as a function of space within the imaging volimie. 

FIG. 5 illustrates one embodiment according to the present invention of a method for 
obtaining image information associated with the dielectric properties of a body from a coil 
array. In step 500, the coil array may be calibrated. Calibration may include various tuning 
and/or adjusting of parameters associated with the coil array as described in greater detail in 
connection with FIG. 7. Generally, calibration is performed in the absence of a load. 
However, in some embodiments, calibration may not be necessary or desirable and may 
therefore be omitted. 

In step 5 10, the coil array may be loaded with the body to be imaged. For example, in 
FIG. 6, body 610 may be placed proximate coil array 600. In step 520, the impedance matrix 
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of the coil array may be measured in the presence of the load at a number of frequencies. By 
operating the coils over a range of frequencies and measuring impedance characteristics, 
information about resonant frequency shifts and other changes in resonant properties may be 
incorporated. For example, an impedance matrix obtained as a frmction of frequency may 
incorporate resonant frequency information that may be used to deteraiine the dielectric 
distribution of an imaging volume. 

The impedance of each coil may be measured using any of various matching networks 
and/or network analyzers as described in more detail m connection with FIG. 9 or using any 
other suitable equipment and/or circuitry. The impedance matrix may be measured by 
operating each coil hi isolation or by operating multiple coils in tandem and measuring induced 
currents, voltages, or other electrical properties in each of the coils. In a tandem operation, 
coils may be operated in pairs, in larger groups, or simultaneously as an entire army. 

In step 530, an unaging volume, including at least a portion of the body, may be 
segmented into a plurality of regions. In general, an image is a set of intensities describing the 
properties of a corresponding region of the imaging volume. Accordingly, a volume 
containing a portion of the body to be imaged may be segmented into a plurality of regions, 
each region having an associated voxel (or pixel) in the resulting image that describes at least 
one intensity of the corresponding region. 

For example, a conductivity image may include voxels representing the value of the 
conductivity of material located in the corresponding region of space. Similarly, a permittivity 
image may be produced, or an image having intensities that incorporate both the conductivity 
and permittivity value of the associated region of the imaging volume into an intensity value. 
Images may also be produced having intensities that represent a broad range of 
electromagnetic properties, including, but not limited to, permittivity, conductivity, 
permeability, electric field, magnetic field and/or any combination of the above. 

The imaging volume may be segmented into regular or irregular regions. For example, 
in FIG. 9, the imagmg volume is segmented into a regular 3D grid of equal sized and 
uniformly distributed cubes. However, one or more regions could be of different size or shape. 
For example, a priori information (e.g., information obtained from another imaging modality, 
separate measurement, knowledge of the body, etc.) may be used to segment the imaging 
volume intelligently to facilitate higher resolution images as discussed in greater detail below. 
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In step 540, an initial trial impedance matrix may be determined. In one embodiment, 
each voxel is assigned an initial trial conductivity and/or permittivity value, that is, each region 
has an associated trial value or values. The trial values may be assigned arbitrarily or they may 
be assigned according to a priori information about the body being imaged. For example, 
regions of suspected homogeneity may be identijBed from other imaging modalities such as an 
MRI image, or estimates of dielectric properties of one or more of the regions (e.g., from 
known characteristics of the body), maybe obtained from a separate source, or obtained from 
computational models approximating the body. 

It should be appreciated that in some enviroiunents (i.e., when certain types of bodies 
are being unaged), permeability may also be expected to vary as a function of space while in 
other bodies permeability may be assigned as the permeability of free space or some other 
constant value. Accordmgly, some embodiments include assignmg trial values of permeability 
to the imaging volume. The trial values assigned to the imaging volume may then be used to 
compute an initial trial impedance matrix. 

For example, the trial conductivity and permittivity values (and in some instances trial 
permeability values) along with the electromagnetic fields computed using Maxwell's 
equations may be employed in the relationship expressed in equation 3 or to compute current 
amplitudes in each of the coils resultmg from known applied voltages in order to compute 
impedance characteristics of the coil array. It should be appreciated tiiat any expression 
describing impedance characteristics as a function of the dielectric properties of the body and 
the electromagnetic fields of the array may be used to determine an initial trial impedance 
matrix. 

In step 550, the initial trial impedance matrix may be compared to the measured 
impedance matrix to determine if the trial values (i,e., the dielectric distribution) assigned to 
the body to form the trial impedance matrix are a satisfactory description of the conductivity 
and/or permittivity distribution of the body as indicated by the measured impedance matrix. 
That is, the trial impedance matrix and the measured impedance matrix may be compared to 
determine if they are sufficiently similar, or that a distance (e.g., a comparison metric) between 
the two is sufficiently small. In one embodiment, the measure of similarity between the trial 
and measured impedance matrix may be a least squares distance, although any measurement 
that indicates similarity between the trial and measured impedance matrices may be suitable. 
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In step 560, a distance between the trial and measured impedance matrix is evaluated. 
If the two matrices are determined to be close enough, that is, if the trial impedance matrix has 
converged to substantially the measured impedance matrix, or the distance between them is 
acceptably small, an image may be formed from the conductivity and/or permittivity values 
used in the formation of the trial impedance matrix. Otherwise, the trial impedance matrix 
may be updated. 

In step 570, when it is determined that the trial impedance matrix and the measured 
impedance matrix are too dissunilar (i-e., a distance between the two matrices is too great) the 
conductivity and/or permittivity values assigned to one or more voxels may be adjusted such 
that an updated trial impedance matrix formed from updated trial values is nearer the measured 
impedance matrix than during the previous iteration. 

Any of various methods for iteratively adjusting the dielectric properties assigned to 
each voxel including, but not Ibnited to, gradient descent, searchmg algorithms such as 
simulated annealing, statistical methods such as expectation maximization, or any other 
optimization method of solving a set of linear or non-linear equations, etc., may be used to 
converge the trial impedance matrix to the measured impedance matrix. Steps 550-570 may be 
repeated xmtil the trial impedance matrix has converged or is determined to be satisfactorily 
close to the measured impedance matrix. In addition, analytical or other non-iterative methods 
may be used for finding a solution. 

In step 565, when the distance between the trial impedance matrix and the measured 
impedance matrix has converged to a local minimimi, an image may be formed from the trial 
values. For example, an image havuig a plurality of voxels may be formed wherem each voxel 
is assigned an intensity value related to the conductivity and/or permittivity value of an 
associated region of the imaging volxmie. 

It should be appreciated that the method described in connection with FIG. 5 depends 
in part on computing an impedance matrix given a set of conductivity and permittivity values. 
Accordingly, results may depend on how accurately the impedance matrix may be computed. 
Applicant has developed methods for calibrating a coil array such that computed impedance 
matrices of the unloaded array agree with measured impedance matrices of the unloaded array. 
A calibrated coil array may give a baseline indication of the accuracy of impedance matrix 
computations. 
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FIG. 7 illustrates one embodiment of calibrating an array of coils according to the 
present invention. For example, the calibration method described below may be used in step 
500 in FIG. 5 to improve the accuracy of impedance matrix computations. Computing an 
impedance matrix may include building a model of a coil array and an imaging volume. 

The term "model" refers generally to a description or representation of one or more 
objects (e.g., a coil array, an imaging volume and its contents, such as a body to be imaged, 
etc.). While a model may emulate a real object (e.g., a model of coil array may have a real coil 
array counterpart), a model is virtual, typically embodied by one or more mathematical 
descriptions. For example, a model may include descriptions of the geometry of the one or 
more objects, parameters tiiat describe the characteristics of the object, etc. Models may be 
stored electronically, for example, on computer memory or as part of a executable program 
stored on computer memory. 

In order that a model emulate at least some properties of the real object, one or more 
functions of the object being modeled may be simulated using the model representation. 
Sunulation refers to computing a function, operation or action such that the model 
representation behaves similarly to its real counterpart. 

Consider the case when the real objects being modeled include a coil array. When one 
or more coils in flie array are operated, electromagnetic fields are generated. This 
electromagnetic environment in turn affects properties of the coil. Simulation of a coil array 
model may include computing the electromagnetic environment of the coil array either in the 
presence or absence of a body acting as a load. Methods and algorithms for performing a 
simulation are often embodied in one or more software programs operating on the model or 
representation of the system being simulated. 

For example, a model of a coil array may be simulated by solving Maxwell's equations 
within an imaging volume defined as part of the model. By solving Maxwell's equations 
during simulation, the electromagnetic environment of an imaging volume may be determined 
such that, together with trial conductivity and permittivity values (during calibration, trial 
values may be set to the conductivity and permittivity values of air), may be used in the 
expression of equation 3 to compute an impedance matrix. 

One method of simulating a model of a coil array (both with and witiiout a load) 
includes solvmg Maxwell's equations according to a Finite Difference Time Domain (FDTD) 
algorithm. In FDTD, a volume of space (e.g., the imaging volume) with or without one or 
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more objects (e.g., a body to be imaged) is partitioned into a lattice or mesh. The 
electromagnetic fields in each region of the mesh are solved for according to an applied set of 
initial conditions (e,g., conductivity and/or permittivity distribution, coil geometry, etc.) and 
boimdary conditions (e.g., the extent of the imaging volume, etc.). 

While FDTD may provide a fast and effective method for simulating the operation of a 
coil array and computing an impedance matrix, any of various other computational methods 
including Chebyshev polynomial expansion, finite element, finite difference frequency domain 
(FDFD) algorithms, any of various other firequency domain computational methods etc., may 
be suitable for simulation and are considered to be within the scope of the invention. 

In step 710, an impedance matrix of the unloaded array may be measured. Various 
methods of measuring the impedance matrix may be suitable to obtain a measurement with 
which to calibrate a model of the unloaded array. In one embodiment, each coil may be 
operated in isolation and the current induced in each of the other coils measured. For example, 
a pulse fi-om an RF power source, or any other electrical stimuli (e.g., a voltage, a current, etc.) 
may be applied at the input of one of the coils and the currents in each of the other coils 
measured. For example, tiie power source may act as a voltage or current source providing a 
broadband electric field pulse. The pulse may be, for example, a "chirp" containing a range of 
frequencies. The range of frequencies may include the nominal resonant frequency of one or 
more of the coils and chosen to include an approximate shift in resonant frequency. This 
process may be repeated at the input of each of the coils in the array. 

In step 720, a model of the unloaded array may be computed. The model may be 
computed at any time and portions may only need to be computed once for a given coil array. 
For example, parameters of the coil array that do not change from one calibration to the next 
(e.g., number and arrangement of coils, various coil properties such as conductivity and 
capacitor values, etc.) may be generated and stored. Parameters such as the size of the imaging 
volume, quantization size of the volume, etc., may need to be input into the model each time 
the imaging environment changes as described in fiirther detail below. 

The army coils and the imaging volume may be partitioned into a lattice, that is, the 
imaging volume may be segmented into a plurality of regions. FIG. 8 illustrates a two coil 
array 800 and a corresponding model of the array 800' shown as a collection of discrete 
squares. The conductive strips 820a-82d may be simulated as a plurality of discrete units (e.g., 
rectangles 820a', 820b', 820c', etc.) having electrical characteristics of a high conductivity 
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material. Likewise, capacitors may be simulated as high permittivity insulating materials 
illustrated by simulated capacitors 830a'-830d'. Each coil may have a simulated port 840'. It 
should be appreciated that other parameters and properties of the array may be defined in the 
model in order for the model to be simulated as desired. 

The imaging volume may be defined as any suitable volume proximate the coil array. 
The imaging volume may then be partitioned into discrete regions. It should be appreciated 
that the imaging volume can be partitioned at any desirable resolution. Other parameters of the 
unaging environment may be defined and added to the model. 

In step 730, an impedance matrix of the unloaded array may be computed by simulating 
the operation of the coil array. For example, FDTD simulations may be performed on the 
model to compute impedance characteristics of each coil in the absence of a load. In step 740, 
the computed impedance matrix may be compared to the measured impedance matrix to 
determine if the model generates values that generally agree with measured values. 

In step 750, parameters of tiie model may be adjusted to compensate for differences in 
the computed and measured impedance matrices. Parameters of the model that may facilitate 
calibrating the model include, but are not limited to, additive correction to the array's 
inductance and/or resistance matrix, adjustments to the nominal value of the matching 
capacitors (e.g., to account for errors in capacitor fabrication and/or stray or parasitic 
capacitance in the matching circuitry), etc., may be adjusted so that the computed impedance 
matrix agrees with the measured impedance matrix as discussed in fiirther detail below. 
Accordingly, a coil array model may be calibrated such that before the array is loaded, 
impedance matrix computations can be made to substantially agree with measured impedance 
matrices of the unloaded array. 

As described above, each region of the body being imaged may have a corresponding 
voxel describing one or more dielectric properties of the associated region. The smaller the 
subdivisions of the imaging volume (i.e., the smaller the regions are chosen) the greater the 
resolution of the resulting image. However, the size of the regions may be limited by the 
resolving power of the coil array. That is, the resolution may depend in part on the coil array's 
ability to ascribe loading eflFects on the properties of the array to individual regions of the 
unaging volume, which depends in part on the number and arrangement of coils in the array. 

The resolution of an RFIM unage may depend in part on the number of correlated or 
independent measurements that may be made of a coil array. For example, the size of an 
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impedance matrix may be related to the resolution of the resulting image. This is, in turn, may 
be related to the nimiber of coils in the array. As the number of coils in an array mcreases 
and/or their size decreases, a size of a region of an imaging volume of which the coil array can 
resolve a conductivity or permittivity value decreases. As such, the size of each region of an 
imaging volume is inversely proportion to the resolution of the image. 

It should be appreciated that an impedance matrix may include both a resistance matrix, 
which encodes primarily conductivity distribution, and a reactance matrix, which encodes both 
conductivity (e.g., via eddy currents) and permittivity (e.g., via capacitive loading). Each of 
these matrices may provide N(N+l)/2 equations, where N is the number of coils. Considering 
redundant information, it is estimated tiiat voxels may be imaged with N coils. However, 
conditioning information, such as information obtained from MRI, other measurements, and/or 
any a priori information about an object to be imaged may greatly improve this number. 

FIG. 9 illustrates one embodunent according to the present invention of an image 
acquisition system for obtaining an RFIM image of a body. Image acquisition system 900 
includes a coil array 910 having a plvirality of RF coils (e.g., coil 910a and 910b). Coil array 
910 is illustrated as having 8 individual coils. However, the arrangement and number of coils 
illustrated in FIG. 9 is merely exemplary and is not limiting, as a coil array may be chosen to 
have any number of coils in any arrangement. For example, a coil array may be arranged such 
that coils are positioned on all sides or a desired number of sides of an object to be imaged 
instead of just underneath a body 950 as shown in FIG. 9. 

Power source 915 may be coupled to coil array 910 to provide power to operate the 
coils. While power source 915 is shown schematically as connected generally to coil array 
910, it should be appreciated that the power source may be connected to the coil array such that 
power (e.g. a voltage or current waveform) may be individually provided to each of the coils in 
the array. That is, each coil may have a separate port with which to receive power. 

Coil array 910 may also include various blocking networks that are capable of 
selectively blocking current flowing aroimd a given coil. Such blocking networks may 
effectively turn one or more coils off, for example, during calibration or during measurement 
of other coils in the array, etc. One or more baluns may also be included to prevent 
unshielded currents from flowing on the cables that connect the array to the other equipment in 
the system. Other circuitry may be included such that each of the coils in the array can be 
selectively operated. 
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Matching circuitry 925 may also be coupled to coil array 910. As with the power 
source, matehmg circuitry is illustrated as generally coupled to the coil array. However, it 
should be appreciated that each coil may have its own matching circuitry such that various 
properties of the coil may be measured by measurement equipment 935 . For example, a 
number of exemplary and suitable matching circuits are illustrated in FIGS. 10a, 10b and 10c. 
Any of various other matching circuits may be used such that properties of the coil array may 
be measured. 

Measurement equipment 935 may be coupled to one or more of the matching circuits to 
measure properties of the coil array such as a voltage, current and/or impedance of mdividual 
coils m the array. Measurement equipment 935 may, for example, include a network analyzer 
having one or more ports connected to the matching circuitry 925 and capable of obtaining 
measurements of one or more properties of the coil array. 

A computer 945 may be coupled to the measurement equipment 935 to receive the 
measurement of one or more properties of the coil array. Computer 945 may be any 
component capable of performmg mathematic computations and/or logic operations. 
Computer 945 may be, for example, one or more microprocessors or digital signal processors. 
Computer 945 may also mclude a computer readable medium such as a memory capable of 
being encoded with instructions, for example, a program configured to perform various 
functions and operations when executed by one or more processors. Computer 945 may be 
included as part of the measurement device 935, for example, a processor included in the 
network analyzer or may be a separate component. 

Computer 945 may be configured to perform computations to facilitate any of 
calibrating the coil array, modeling the coil array, computing impedances matrices, forming an 
image of dielectric properties of a body, etc. Computer 945 may also be coupled to a display 
955 capable of rendering RFIM images acquired from measurements of the coil array. 

FIG. 1 1 illustrates one embodiment according to the present invention of forming an 
RFIM image of a body using an image acquisition system, for example, image acquisition 
system 900 illustrated in FIG. 9. In step 1 100, the coil array may be calibrated. During 
calibration, power source 915 in FIG. 9, for example, may apply a waveform at a desired 
frequency (e.g., an RF pulse) to the coil array such that a current is generated in one of the 
coils. As discussed above, an impedance matrix may be formed from S-parameters measured 
by, for example, a network analyzer. 
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S-parameters describe the transmission and reflection of traveling waves, that is, they 
may represent the reflection and transmission coefficients between incident and reflection 
waves and are often measured as a function of frequency. The S-parameters may be used to 
describe the behavior and characteristics of circuit or network (e.g., a resonant coil). There 
exists a one to one correspondence between S-parameters and unpedance and the two 
measurements can easily be converted into one another, for example, by lS'=(2^Z^^;fZ+Za;'^ 
where Zq is the match impedance of the RF source (e.g., 50 or 75 ohm BNC impedance). 

Accordingly, a measured impedance matrix may be obtamed by successively operating 
each coil at a range of different frequencies and measuring the S-parameters of the coil array. 
That is, power source 915, may supply power to each coil over a range of desu^d frequencies. 
Typically, the range of frequencies includes a nominal resonant frequency of the coil being 
operated. For example, a coil may be designed to resonate at a given frequency. As discussed 
above, when a coil is loaded, the resonant frequency may shift in response to the 
electromagnetic coupling between the coil and the load. 

When a coil is operated at its resonant frequency, a maximum amount of energy is 
being coupled into the coil and a minimum amount of energy is being reflected back to the 
power source. This resonant condition manifests itself in a local minimum in the magnitude of 
the so-called Sn parameter. Accordingly, if a given coil is operated at a range of frequencies 
including its nominal resonant frequency and the Su parameter measured at each desired 
frequency in the range, a plot of the Su parameters will show a dip in magnitude at the 
resonant frequency. When a coil is loaded and the measurements repeated, the dip in 
magnitude will have shifted due to the shift in resonant frequency of the coil. As such, 
resonant frequency shifts may be encoded in S-parameter measurements. In addition, changes 
in depth and breadtii (in frequency) of this dip encode further information about the resonance. 

In addition to the resonant frequency shift mformation encoded in the S-parameters, 
resistive, capacitive and inductive coupling between the coils may also be available from 
measurements of the S-parameters. For example, when a given coil is operated, some of the 
power is dissipated in the other coils. This power can be measured in the S-parameters of the 
other coils. The notation Sij will be used to indicate the S-parameter measured in the f' coil in 
response to operating the i'* coil. 

Accordingly, an unpedance matrix may be obtained as Sj/co), where © are the 
frequencies over which the coil is operated and S-parameters of the array are measured. This 
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measured impedance matrix may incorporate information about resonant fiequency 
perturbation, resistive, capacitive and inductive coupling between the various coils in ihe array. 
This information may be employed to determined tiie dielectric properties of the loading body. 

In step 1 120, a model 1 125 of the coil array may be constructed. In particular, a model 
may be defined such that the electromagnetic environment of the coil array resulting from 
operating an arbitrary number and arrangement of coils may be simulated. The geometry, 
properties and characteristics of tiie coil array may be described (e.g., as a matiiematical 
representation) so that the operaticm of the coil array may be simulated. For example, an 
operator may program a model of the coil array by describing various operating parameters, 
values of vutual components (e.g., capacitor values, eto.), geometry of the array (e.g., number 
of coils and spatial relationship between coils), etc., and store the model on a computer 
readable medium, such as a memory stored, for example, on computer 945 in FIG. 9. 

In addition, an imaging volume wherein a dielectric distribvition will be computed may 
be defined. For example, in FIG. 9, imaging volume 905 may be defined and included in tiie 
model. For example, imaging volume 905 may be chosen to at least enclose a portion of a 
body tiiat is desired to be imaged. The imaging volume may be partitioned into a plurality of 
regions (e.g., such as region 905a, 905b, 905c, etc). 

It should be appreciated that the imaging volume may be partitioned in any number of 
ways and the volume illustrated is merely Kcemplaiy. For example, the size of each region 
may be chosen to provide a desired resolution of the resultmg unage. In addition, the regions 
may be chosen to be of different size and shape as suits a particular implementation (and in 
consideration of the body being imaged). The imaging volume may be chosen to enclose 
certMn portions of the object, may be chosen to enclose the entire body, or chosen to enclose 
the body and the coil array and/or any peripheral space around tiie coil array that facilitates 
simulation of tiie coil array. The constructed model 1 125 may be stored and later accessed 
during coil array simulation. 

In step 1 130, a computed impedance matrix may be generated. During the computation 
of a computed impedance matrut for calibration, each of the plurality of regions in the imaging 
volume may be assigned values corresponding to the dielectric properties of air to simulate an 
unloaded coil array. For example, conductivity values ao and eo may be assigned to the 
plurality of regions, where ao and eo are the conductivity and permittivity values of the empty 
imaging volume (e.g., ah). 
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The model of the unloaded coil array together with the assigned dielectric properties of 
the empty imaging volume may be simulated, for example, by performmg an FDTD simulation 
of the model. It should be appreciated that the model, the imaging volume, and the FDTD 
algorithm (in addition to other functionality) may be stored on a memory of computer 945, for 
example, as one or more software programs. 

The simulation may compute an impedance matrfac by determinmg the unloaded EM 
envu-onment and, together with the assigned dielectric values, employ the relationship 
expressed in equation 3, In particular, smce the conductivity and permittivity values have been 
assigned as a discrete function of space (i.e., the dielectric properties assigned to each region of 
the imaging volume), and the electric and magnetic fields generated by the coil array at various 
input frequencies may be calculated throughout the imaging volume (e.g., by solving 
Maxwell's equations during shnulation), the impedance values of the coils may be computed. 

Radiation losses may also be included in the simulation (even though the effects may 
be relatively small). In one embodmient, computation of radiation losses may be achieved in 
an FDTD simulation by choosing the appropriate boundary conditions. For example, the 
appropriate 'out-gomg wave' boundary conditions may be applied at the surface extent of the 
unagmg volmne. Accordingly, a computed unpedance matrix and a measured impedance 
matrix may be obtained for the unloaded coil array. 

In step 1 140, this information may be used to adjust parameters of the model such that 
computation of impedance characteristics agree with measurement of unpedance 
chamcteristics. Stated differentiy, the model may be adjusted so that shnulation of the coil 
array comports with the actual operation of Ihe coil army. 

For example, additive correction may be made to the array's resistive, capacitive and/or 
inductive matrix. In some embodhnents, this correction may mitigate errors that may arise due 
to the discrete nature of the simulation of the model. For example, in an FDTD simulation, the 
partitioned regions of the imagmg volumes may not take into account, for example, the 
contribution to an array's impedance at distances very near the coils (i.e., at distances smaller 
than the quantization of the imaging volume). In addition, defects or imperfections in the 
fabrication of the coils may need to be compensated for by adjustmg the parameters of the 
model. Errors m capacitor fabrication and stray capacitance m, for example, tiie matching 
cu-cuitry may be accounted for as well. Otiier adjustments may be made as well to compensate 
for differences in the model and the real coil array being represented by the model. 
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When a network analyzer is calibrated, the lengfli of a cable connecting the network 
analyzer to the coil army may be taken into consideration. However, small deviations from the 
assumed length of the cable (e.g., due to small additional circuitry in and near the connector, 
such as a BNC connector) may add unaccounted for phase to the measured S-parameters, 
which may cause the computed S-parameters and measured S-parameters to disagree (i.e., 
phase shifted from one another as a function of the cable length deviation at each coil port). 
This phase may be accounted for m the model by shifting the phase of the measured or 
computed S-parameters accordingly. For example, phase correction may take the form: 

S,{a>) >e''^^-y^^'^%{c>)^ (6) 

where q)i(a,) is the phase acquired at coil port / at frequency co (such phases increase 
linearly with frequency). The various adjusted parameters of the model may then be used to 
update model 1125. 

In step 1 150, the coil array may be loaded with a body to be imaged. For, example, 
body 950 in FIG. 9 may be placed proxhnate the coil array 900 such that the loading effect of 
the body can be measured, for example, m the loaded S-parameters of the coil array. 

In step 1 160, a measured impedance matrfac for the loaded array may be measured. 
Measuring the impedance matrix witii the body proximate the coil array may be performed in a 
manner similar to that described in connection with measuring the impedance matrix of the 
unloaded array. For example, in the system of FIG. 9, each coil may be operated successively 
over a range of frequencies. At each frequency, the S-parameters for the coil array may be 
obtained via matching circuitry 925 and measurement device 935. 

In step 1 170, the model may be updated according to an initial approximation of the 
dielectric distribution of the imaging volume with the addition of the loading body. For 
example, the imagmg volume may be assigned initial trial conductivity and/or permittivity 
values. The trial values may be assigned arbitrarily or may be assigned according to a priori 
information about the imaging body. Any a priori information available that may be employed 
to condition the simulation may facilitate faster convergence times and may reduce 
opportunities for the sunulation to converge to an undesirable local mmunum, and may result 
in the ability to obtain higher resolution images. 
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In step 1 180, operation of the coil array may be simulated according to the model 1 125 
to produce a trial impedance matrix. Simulation may be performed substantially as described 
above for calibration using the updated dielectric properties of the imaging volume. In step 
1 190, the trial impedance matrix may be compared with the measured impedance matrix to 
determine a distance between the two matrices. This comparison may provide an indication of 
how closely the trial values approximate the actual distribution of the dielectric properties of 
the imaging volume. 

If the distance is determined to be too great, that is, if the distance suggests that the 
distribution of dielectric properties is not a close enough approximation of tihie true dielectric 
distribution of the body, tiie trial values of conductivity and/or permittivity may be updated. 
For example, the trial values may be updated according to a method of gradient descent. 
According to one method, the system of equations expressed as, 

Z^^'(^,^,„^„) = Z--(^), (7) 

may be solved, where (ij) indicates the coil pah: index and the index n enumerates the 
regions into which the imaging volume is divided. Solvmg the system of equations expressed 
in equation 7 may be achieved in a number of ways. In one embodiment, the Port optimization 
library developed by AT&T is used, although various other optimization schemes may be 
employed. 

As is often the case with optimization schemes, the inquiry of whether a problem is 
well-posed of often addressed, for example, whether the system of equations expressed in 
equation 7 is well-posed. This inquiry typically involves asking whether a solution exists, and 
if so, if the solution is unique. In the RFDM framework as described above, these questions 
may be answered by considering the so-called "electrical prospection" problem. 

Consider an imaging volimie V. By specifying a tangential component Et of the 
electric field on the boundary of V, Maxwell's equations may be used to compute the electric 
and magnetic fields throughout the entirety of V. In particular, the normal component Bn of 
the magnetic field on the boundary of V may be computed. This constitutes a map of the form: 



(8) 
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for each point on the boundary of V. In RPIM, Et may be controlled by means of the 
coils. For instance, consider an imaging volume whose surface is uniformly tiled by a very 
large number of rectangular coils. By applying known voltages at the mput port of each coil, 
tangential electric fields of known magnitude may be generated. The normal component of the 
magnetic field B^may then be obtained by measurmg the current induced in each of the coils 
via Faraday's law of induction. By using an array that consists of an infinite number of 
infinitely small coil elements and that fully encloses the volume V,.the fiinction F in equation 8 
may be fiiUy determined. It has been shown that the measurement of this function is sufficient 
to uniquely determine local conductivity and permittivity inside V. 

Smce, in practice, only a finite number of coils are available with which to test the 
function F appearmg in equation 8, there may be pathological instances where the inverse 
problem is ill-posed. However, such pathological cases are common to all techniques of 
medical imaghig and remote sensing in general. In practice, these ill-posed cases result in 
image artifacts that a tramed practitioner can recognize and correct and/or that can be avoided 
by appropriate array design and/or conditioning the reconstruction of the image. 

Steps 1 170-1190 may be repeated until it is determined that the trial dielectric 
distribution is a satisfactory approximation of the dielectric distribution within the imaging 
volume, and hence, the dielectric properties of at least a portion of the body bemg imaged. 
That is, the distance between a final trial impedance matrix and the measured impedance 
matrix is suflSciently small (or it has been determmed that the algorithm has converged or 
both). 

In step 1 195, the conductivity and/or permittivity values assigned to the image volume 
to compute the final trial impedance matrix may be used to form an RFIM image. For 
example, each voxel in the image may have an intensity related to a conductivity and/or 
permittivity value assigned to an associated region of the imaging volume. A conductivity and 
a permittivity hnage may be formed, or a single image may be formed deriving intensity values 
fi-om some combination of the conductivity and permittivity values. The resulting RFIM 
image may be stored in memory or displayed, for example, on display 955 as shown in FIG. 9. 

It should be appreciated that a coil array may comprise any number of coils in any 
arrangement. For example, FIGS. 12A and 12B illustrate one embodiment according to the 
present mvention of a coil array arrangement having a plurality of coils arranged substantially 
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in a rectangle configuration. FIG. 12B iUustrates the same array as in FIG. 12A, however, all 
of the coils in the array are shown using dotted lines for coils on flie hidden feces. 

In coil array 1200, opposing feces of the rectangular array have coils that may be 
positioned essentially orthogonal to one another. In particular, generaUy rectangular coils on 
opposing feces may have tiie axis along Aeir larger dimension positioned essentially 
perpendicular to one another. In other embodiments, coils on opposing feces are arranged 
essentially paraUel to one another. Other non-parallel or non-orthogonal arrangements may be 
suitable to form a generally closed coil array. Coil array 1200 may, for example, define the 
outer boundary of an imaging volume or the imaging volume may be chosen to enclose the coil 
array. 

FIG. 13 illustrates another embodiment of a coil array according to tiie present 
invention. Coil array 1300 may be substantially similar to coil array 1200. However, coils on 
two of the feces may be removed such that the coil array has an openmg. Coil array 13 may be 
well suited for imaging the torso of a patient or portions of an object that may not fit entirely 
within a closed array. 

Coils in an array may be arranged in any ofter arrangement such as in a generally 
cylindrical configuration around a body to be imaged or any other layout wherein loading 
effects of a body may be measured. By increasing Ae number of coils in the array, the 
resolution of an RFIM image may be increased. In addition to the number of coils, the 
arrangement of coils may affect the resolution of image, in part, due to conditioning of the 
imaging environment. 

hi one embodiment according to the present invention, a coil array layout may be 
chosen in order to increase the amount of information that may be encoded about a loadmg 
body. In some circumstances, information about a body to be imaged may be available. For 
example, m a human body, approximate size and shape of the body may be available. In 
addition, nominal values for conductivity and permittivity may be available. This information 
may be used to evaluate the conditioning of a particular array layout. 

In some embodiments discussed above, a system of equations may be solved in order 
to determine the conductivity and permittivity distribution within an imaging volume. For 
example, the system of equations expressed m equation 7 may be solved. This system of 
equations is generally a non-linear system of equations. However, the system of equations 
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may be approximated by a linear system of equations such that linear techniques may be used 
to determine the conditioning of the system of equations. 

For example, ameasure of how well conditioned a system of linear equations is may be 
determined from the condition number. Consider a system of linear equations expressed as Ax 
= 6, where it is desked to solve for x in terms of b with a coefficient matrix^. When the 
coefficient matrix A has a singular value of zero, x cannot generally be solved for because A'^ 
does not exist. In the numerical framework, singular values close to zero may prevent solving 
for X in a shnilar manner. Arrangements in which flie array layout is not well conditioned lead 
to these degenerate or near degenerate conditions. 

In order to avoid such coil layouts, an array may be arranged such that the coefficient 
matrix (i,e., a system of equations) obtained from, for example, simulating the coil array with 
approximations to the make-up of the body (e.g., size, shape and dielectric distribution) does 
not have zero or near zero singular values. One measure may include computing a ratio of a 
largest singular value to a smallest singular value. When the smallest singular value is small 
compared to the largest smgular value (i.e., when the ratio is large) a coil layout may not be 
well conditioned and may need to be adjusted. Accordingly, different coil layouts may be 
arranged such that layout is suitably conditioned. 

In one embodiment, the Jacobian of the non-lmear system of equations is used as a 
Imear approximation of the non-linear system of equations. The singular values of the 
Jacobian matrix may then be determined analytically or by any of various methods such as the 
Jacobi method, Gauss-Seidel, Singular Value Decomposition (SVD), or by any other method 
suitable for obtaining the singular values of the Jacobian matrix. The singular values may be 
compared m any number of ways to determine whether the layout is well conditioned. This 
process may be repeated until a layout with suitable conditioning is obtained. 

In another embodiment according to the present invention, a single coil may be used to 
image a loadmg body. For example, the spatial relationship between the single coil and the 
body may be varied and measurements of the loading effect taken at each placement. It should 
be appreciated from the foregoing that the Sn parameter measures the reflected energy in a coil 
resulting from operating the same coil. Accordingly, this information may be available absent 
interactions with any other resonant coils to determine the dielectric distribution of the body. 

Furthermore, one or more coils may be added to increase the amount of information 
available (e.g., resistive, capacitive, and/or inductive coupling between coils) for determining 
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the dielectric properties of the body. These additional one or more coils may be varied relative 
to the body such that multiple measurements may be taken. By varying the spatial relationship 
between one or more coils and the loading body m relatively small increments, higher 
resolution images may be achieved without requiring relatively large numbers of coils to be 
employed in an array. 

Employing relsonant coils (as opposed to inductive coils such as solenoid coils) 
facilitates using the sensitive loading properties of resonant coils to determine the dielectric 
distribution of a volume of space. In addition, by considering a fuller electromagnetic 
environment of the volume (e.g., resistive, capacitive and/or inductive coupling), more 
mformation is available tiiat may resuft m higher quality, better resolution raiages. 

It should be appreciated that various embodiments according to the present invention 
result in relatively inexpensive and portable equipment. While MRI coils may be used to 
practice aspects of the present invention, simpler and less expensive resonant coils may also be 
appropriate. For example, various equipment symmetries, tuning and detuning circuitry, 
gradient coils, etc., needed for MRI are not required in RFBM. 

In addition, e^qjensive and bulky magnets required for MRI imaging may be 
elimmated. Accordingly, it may be relatively inexpensive to produce an RFIM imaging 
system. Moreover, the equipment has a relatively small footprint and may be capable of being 
transported from one location to another, making it suitable for mobile and/or emergency 
situations. 

It should be appreciated that various aspects of the present invention may be may be 
used alone, m combination, or in a variety of arrangements not specifically discussed in the 
embodiments described in the foregoing and is therefore not limited in its application to the 
details and arrangement of components set forth in the foregoing description or illustrated in 
the drawings. The mvention is capable of other embodiments and of being practiced or of 
being carried out in various ways. In particular, various aspects of the present invention may 
be practiced with any number of coil types and arrangements. For example, generally planar 
coils, birdcage coils, surfece and volume coils may be used alone or in any combination with 
the any of the various imaging techniques described herein. 

In addition, various aspects of the invention described in one embodiment may be used 
in combination with other embodiments and is not limited by the arrangements and 
combinations of features specifically described herein. Various alterations, modifications, and 
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improvements will readily occur to those skilled in the art. Such alterations, modifications, 
and improvements are intended to be part of this disclosure, and are intended to be within the 
spirit and scope of the invention. Accordingly, the foregoing description and drawings are by 
way of example only. 

Also, the phraseology and terminology used herein is for the purpose of description and 
should not be regarded as limiting. The use of "including," "comprising," or "having," 
"containing", "involving", and variations thereof herein, is meant to encompass the items 
listed thereafter and equivalents thereof as well as additional items. 

What is claimed is: 



